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Abstract 


The  problem  of  designing  a  feedback  compensator  to 
minimize  a  weighted  L*  norm  of  the  sensitivity  function 
of  a  M1M0  linear  time  invariant  system  is  considered. 
Tne  problem  is  solved  by  establishing  its  equivalence 
to  the  different  but  related  problem  of  multivariable 
zeroeth  order  optimal  Hankel  approximation  solved  re¬ 
cently  by  Rung  and  Lin. 

1.  Introduction 

In  this  paper  the  oroblem  of  designing  a  feedback 
compensator  to  minimize  the  sensitivity  function  of  a 
M3.V0  linear  time  invariant  system  is  considered.  Sen¬ 
sitivity  is  measured  by  a  weighted  L*  norm.  The  use  of 
a  weighted  L“  norm  to  measure  sensitivity  was  first 
proposed  by  Zames  [ij.  Za'-'es  has  argued  that  a 
weighted  L“  nonn  arises  naturally  as  the  optimization 
criterion  in  sensitivity  minimization  problems  involv¬ 
ing  disturbances  with  variable  but  founded  power  spec¬ 
tra.  In  contrast,  the  quadratic  norm  used  in  the 
Wiener-Hopf  approach  is  a  meaningful  criterion  only 
when  the  disturbances  have  a  fixed  powe-  spectrum  [2]. 
For  S1S0  systems  the  problem  of  sensitivity  minimiza¬ 
tion  in  the  L"  setting  has  been  solve:  ty  Zames  and 
Francis  [3-<].  Safonov  end  Cnen  [s]  provide  a  solution 
for  the  M1M0  case  when  the  sensitivity  is  constrained 
to  be  decoupled  (i.e.,  diagonal).  Francis,  Helton  and 
Zames  [6]  and  Chang  and  Pearson  [7]  have  solved  the 
•  MIMO  problem  without  decoupling  constraints.  Here  we 
.present  a  different  solution  for  the  MIMO  case  which 
) establishes  important  links  between  L*  sensitivity  op- 
vtimization  and  Hankel  norm  optimal  approximation  theory. 


■  Fig.  1 

I  For  the  feedback  configuration  of  Fig.  1,  the  sen- 
Jsitivity  function  is  S  ■  (I* PC)'*  which  is  the  trans- 
•  fer  function  from  the  reference  signal  r  to  the  error 
Isignal  e.  The  weighted  sensitivity  function  we  consid¬ 
er  is  T  +  S  Wr*T  where  Wx  and  Wr  are  appropriately 
chosen  weighting  matrices.  The  optimizction  problem  is 
solved  in  essentially  three  stages:  1)  The  class  of 
sensitivity  functions  for  which  the  closed  loop  system 
is  stable  is  parameterized  by  the  fracfonal  represen¬ 
tation  approach  [6].  2)  The  problem  is  translated  in¬ 
to  one  of  optimal  tT  Interpolation.  3)  A  solution  to 
the  interpolation  problem  is  given  by  establishing  its 
equivalence  to  the  different  but  related  problem  of 
zeroeth  order  optimal  Hankel  norm  model  approximation 
solved  by  Kung  and  Lin  [9,10],  who  extended  the  earlier 
S1S0  results  of  Adamjn,  Arov  and  Krein  [11]. 

Notation  Is  Introduced  in  section  2.  In  section  3 
realizable  sensitivity  functions  are  defined  and  para¬ 
meterized.  In  section  <  the  sensitivity  optimization 
problem  is  shown  to  be  equivalent  to  an  optimal  Inter¬ 
polation  problem.  Its  links  with  the  optimal  Hankel 
norm  approximation  problem  are  established  in  section 
S.  Concluding  remarks  are  presented  in  section  6. 

2.  Notation  and  Preliminaries 

Me  let  H**  denote  the  space  of  matrices  whose  ele- 
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ments  are  functions  analytic  and  bounded  in  the  r.h.p. 
and  let  L*  denote  the  space  of  matrices  whose  elements 
are  functions  bounded  on  the  jw  -  axis.  The  space  H“ 
is  a  subspace  of  L*.  The  norm  of  6(s)  €  L*  is  defined 

as 

S<$.  *  sup  o(G(jw))  (2.1) 

u 

L.  * 

where  F(A)  denotes  the  largest  singular  value  of  A.  A 
matrix  function  is  scdbU  if  its  elements  are  analytic 
in  Re(s)  a  0.  The  elements  of  H“  are  stable  and  proper. 
A  matrix  A(s)  is  all-paaa  if 

A*( jw)  A(ju>)  -  I  for  all  u  (2.2) 

where  A*(s)  denotes  A^(-s).  A(s)  is  inner  if  it  is  all 
pass  end  stable.  B(s)  is  min.  phase  If  it  has  no  r.h.p. 
zeros  and  outer  if  it  is  min-phase,  stable  and  proper. 

Two  elements  of  a  ring  are  right  (left)  coprime  if 
their  only  common  right  (left)  factors  in  the  ring  have 
Inverses  in  the  ring.  Me  will  apply  the  concept  of  co- 
primeness  to  the  rings  of  rational,  stable  and  rational, 
stable,  proper  matrices. 

The  plant  and  the  compensator  transfer  functions 
are  denoted  by  P  and  C  respectively.  P  is  assumed  to 
be  an  nxn  real,  rational,  proper  matrix  having  m  i  0 

distinct  r.h.p.  zeros  zi,  Z2 . Zip.  and  no  poles  or 

zeros  on  the  jw  -  axis.  The  sensitivity  function  is 
S  «  (1  +  PC)*1  and  the  weighted  sensitivity  function  is 
T  *  Wx  $  WV*1.  The  r.etrices  Wx  and  Ur  reflect  relative 
frequency  weightings  on  the  sensitivity  function  and 
the  reference  signals  in  the  sense  that 

IK  S  W  “’a.  «  sup  SWX  5  «!l2 
°/x€L2  ^x,l2 

where  J|-J|2  denotes  the  L*  norm.  We  assume  Wx  and  W,-*1 
to  be  rational,  stable,  proper  and  min.  phase  and  of 
unit  norm.  The  stable,  min.  phase  condition  is  not  a 
restriction  since  if,  sey,  H/(s)  is  non-min.  phase  t^en 
we  can  take  in  its  place  the  min.  phase  spectral  factor 
of  WxT(-s)  Wx(s)  which  gives  rise  to  the  same  frequency 
weighting  as  Wx(s)  [3,  12]. 

3.  Realizable  Sensitivity  Functions 

In  this  section  we  characterize  those  sensitivity  func¬ 
tions  which  correspond  to  a  stable  closed-loop  (c.l.) 
system. 

Oef:  A  rational1  sensitivity  function  S  Is  realizable 
T7T  S  *  (I  +  PC)*1  for  some  C  which  stabilizes  the 
closed-loop  system. 

Lemma  1:  There  exist  rational,  stable,  matrices  Br, 
Dr^  Nx,  Ox,  Ui  and  Vx.with  Br,  Ox  inner  and  U',  V* 
proper  such  that  *  * 


R  U  »  •  i-  B 

r  *  a  r 
Kr  P  -  01*1  Nj 

«XUX  4  Vi  "  1 


•*(s)  4  (s  -  z,) . 


Proof:  Since  P  Is  a  rational,  proper  matrix,  it  admits 
of  right  and  left  fractional  representations 

^All  matrices  used  in  this  section  are  rational. 


p  “  Nrl  °r1  *  °i1  Nil  (3-^ 

where  Nri,  Dri  are  right  coprime  and  N/i ,  0/i  are  left 
coprime  in  the  ring  of  rational,  stable,  proper  ma¬ 
trices  [13].  W/  Nri  is  rational,  stable,  proper  and 
can  be  factored  as  We  Nri  «  Aj  Ao,  where  A.j  is  inner 

and  An  is  outer  [12].  Define  a*(s)  A  li  (s  -  z,)  and 

1*1 

n  .  a*{$)  m  -1  Tv  _  _  k  v - _ _  v  £  a.  Lv,.«  , Vv 


Br  l  *44 
r  =  a  ( s ) 


The  r.h.p.  zeros  of  A*,  being  the 


r.h.p.  zeros  of  Nri  and  in  turn  the  r.h.p.  zeros  of  P, 
are  contained  in  a*(s).  Hence  Br  is  stable.  Further, 

—  and  A,  being  all-pass  8r  is  all-pass.  Hence,  8r  is 

°  '  a*  -1 

inner.  We  can  now  write  W/  Nri  «  —  Br  Ap.  Aq  being 

outer  we  can  factor  Dri  *  Dr  Aq,  Dr  stable.  Then 

Br  P  *  V  “r*1- 

■ 

Similary  factor  D/i  Wr  *  A^-D/,  where  AqI  is 
outer  and  0/  is  inner.  Again,  Aq]  being  outer  we  can 
factor  Njjj  =  Agi-N/,  Njj  stable,  then  Wr  P  *  D/*'N/. 

Since  N/i,  D/i  are  left  coprime  in  the  ring  of 
rational,  stable,  proper  matrices,  there  exist  rational, 
stable,  proper  U/i ,  V/i  which  satisfy  the  Bezout  iden¬ 
tity 

Nil  UX1  +  °A1  VZ1  5  1  (3  2) 

[13].  If  we  substitute  N/i  *  Agi  N/  and  D^i  *  AqI  0/  Hr 
then  Ajji  N/  U/i  +  AqI  D/  Wr  V/i  =  I  which  is  equivalent 
to  N /  U/i  AqI  ♦  D,  Hr  V/i  AqI  *  1.  Now  U/  £  U/I-AqI, 
stable,  proper  and  V/  ^  Wr  V/i  Aqi  ,  stable,  satisfy 

",  V  »,»,•'  13-3> 

Eq.  (3.3)  implies  that  D/  V/  is  proper.  Since  0/  is 
inner  this  implies  that  V/  is  proper.  Q.E.D. 

Leima  2:  S  is  realizable  iff 

°  L’  V  -  ^  B/1  X)D,  (3.4) 

X  o  r  X 


T  i  Hx  S  Ur* 


(UxHr 


alent  to 
T  > 


w  S  w 
i  r 


¥  Br*1  X)Di 


subject  to  the  constraint  that  the  closed-loop  system 
of  Fig.  1  is  stable.  o 

As  a  consequence  of  Lemna  2,  problem  1  Involves 
minimization  over  the  set  of  weighted  sensitivity  func¬ 
tions  parameterized  as 

T  -  (Wx  W,.*1  Vx  -  ^  B/1  X)0X  (4.2) 

X  rational,  stable.  If  we  take  X  *  0  then  ^W/W^VD/ 
is  proper  and  has  finite  norm.  Consequently,  the  minl- 
mun  attained  in  probleml  is  finite.  Also,  jjTJj,  «  -  if 
T  is  improper.  Hence  we  can  restrict  ourselves  to  mini¬ 
mizing  over  proper  T's  only.  We  can  also  relax  the 
constraint  that  T  be  rational  since  the  optimum  T  turns 
out  to  be  rational  [4,5].  With  these  modifications  and 
recalling  that  stable,  proper  T's  are  parameterized  by 
X  €  H-  we  have  an  equivalent  problem. 

Problem  2; 

Min  |JT(X);U  (4.3) 

X  €  H" 


T(x )  *  (WxWr* 


4  B/1  X  )D, 


If  we  define 


*  k  \  \  Wr 
ilK-yll 


Y  =  Br  T  Dx‘ 


for  some  rational,  stable  X. 

Proof:  Note  that  (3.1)  is  a  coprime  fractional  repre¬ 
sentation  of  P  and  (3.2)  is  the  corresponding  left 
Bezout  identity  in  the  ring  of  rational,  stable  ma¬ 
trices  also.  Hence  from  [8],  S  is  realizable  iff 

S  *  (NrlZ  +  V/l)D/l  (3'5) 

for  some  rational,  stable  Z.  If  we  substitute  for  Nri , 
V/i  and  D/i  from  the  proof  of  Lenina  1,  (3.5)  is  equiv¬ 
alent  to 

T  *  WXS  *  (WxWr-1  Vx  -  £  B/1  X)DX 

for  X  *  -  B0Z  AqI.  B0  and  A^  being  rational,  outer, 

X  is  rational,  stable  iff  Z  is  rational,  stable.  Hence, 
5  is  given  by  (3.5)  for  some  rational,  stable  Z  iff  T 
is  given  by  (3.4)  for  some  rational,  stable  X  .  Q.E.D. 

Eq.  (3.4)  gives  a  parametrization  of  rational, 
stable  realizable  weighted  sensitivity  functions  T  in 
terms  of  rational,  stable  X.  In  the  following  section 
we  will  also  consider  T's  which  are  given  by  (3.4)  but 
which  are  rational,  stable  and  proper  or  only  stable 
and  proper.  The  former  are  parameterized  by  rational, 
stable,  proper  X  and  the  latter  by  stable,  proper  X, 
1.e..byX€H».  <1.E-D- 

4.  The  Interpolation  Problem 

In  this  section  we  show  that  the  characterization 
(3.4)  of  realizable  sensitivity  functions  translates 
the  sensitivity  minimization  problem  into  one  of  opti¬ 
mal  If*- interpolation.  Our  objective  is  to  solve 

Problem  1:  Find  a  rational  compensator  C(s)  to  mini¬ 
mize 

HWAd  +  PC)'1  V1!!.  (4.1) 


Since  Br,  D^  are  inner 

l|Y||.  •  JTJ.  (4.7) 

and  problem  2  is  equivalent  to 
Problem  3: 

Min  !|K  -  £  x;|.  (4.8) 

X€H*  “  □ 

If  Yb  achieves  the  minimum  in  (4.8),  then  the  op¬ 
timal  sensitivity  is 

S°  *  (B  W.)'1  Y°  (D.  Wj  (4.9) 

I  i  At  i 

and  the  optimizing  compensator  is 

r°  -  P'^S0'1  '*  1r" 


(4.10) 


Remark:  Problem  3  is  solved  by  Sarascn  [15]  in  t*e 
context  of  ^-interpolation.  Francis  and  Zames  [4] 
show  the  sensitivity  minimization  problem  to  be  equiva¬ 
lent  to  problem  3  for  the  SISO  case  and  Chang  and 
Pearson  [7]  for  the  MIMO  case.  In  both  cases  Sarason's 

results  are  finally  used  to  obtain  the  minimal  sensi¬ 

tivity.  Here  we  take  a  different  approach.  We  convert 
problem  3  into  a  Hankel-norm  approximation  problem  and 
thus  establish  an  important  link  between  1“  sensitivity 
minimization  and  optimal  Hankel  approximation. 

The  following  Lemma  shows  that  the  requirement 

that  Y  be  of  the  form  Y  *  K  -  ~  X  for  some  X  €  if1  is 
In  fact  an  interpolation  constraint  on  Y. 

Lemna  3:  (a)  V  *  K  -  ~  X  for  some  X  g  H*  iff  Y (z , )  * 
K(zi ) ,  1  *  1  *  m. 

(b)  Y°  *  Mfn{3Vj|-|Y  -  K  -  ¥*•  X  6  H"}  (4  11} 

*  M1n{J| Y’l.l  Y  -  A  -  y*,  X  6  H-) 

for  only  A  satisfying  A (zi)  «  K(Zf),  1  *  1  <  m. 

Proof:  (a)  If  Y  *  K  -  yX  for  some  X  €  H"  then  for 
Itis  m,Y(zt)  •  K(zi)  sine;  a*( z ^ )  *  0.  Conversely  , 
If  Y(zi)  «  <C(zi)  for  1  <  i  <  a  then  yr(Y-K)  is  stable 
and  proper,  i.e.,  js(Y-K)  €  iT.  Hence  Y  *  K-yX  for 


soar  X  €  IT. 

(b)  From  part  (*),  Y  ■  K  -  for  some  X  €  IT  Iff 
Y(z< )  *  K(z,),  1  *  1  *  m.  Similarly  Y  ■  A  -  ^pX,  for 
some  X,  €  H-  iff  Y(z,)  ■  A(z,).  Is  1  <  a.  But  since 
A(zi )  •  K(zi),  lsisi,  Y  ■  K  -  ~X  for  some  X€»T 
Iff  Y  *  A  -  YX1  for  iame  X1  *  Eq.  (4.11)  then 
follows. 

As  a  consequence  of  Lenina  3,  problem  3  Is  essen¬ 
tially  an  optimal  IC-interpolation  problem,  l.e. ,  one 
of  finding  a  Y  €  IT  of  minimum  norm  which  satisfies 
the  r.h.p.  Interpolation  constraints  Y(zy)  *  K(zy), 
l*l*i. 

Me  next  present  the  Interpolation  problem  In  a 
form  which  leads  to  a  link  with  a  Hankel  approximation 
problem. 

Theorem  1;  Problem  3  has  the  same  solution  Y°  as 
Problem  4: 

Min{:|Y||.|Y  «  ^<N*  -  X),  X  €  H")  (4.12) 


m 

H*(s)  ti  l 
*  1  «1 


Theorem  2:  if  Ha®(s)  Is  an  optlMl  zeroeth  oroer  Han- 
kel  approxlmant  of  H(s)  and  Eo  the  corresponding  error, 
then 

(a)  *in  ||£(H*-  JOB.  •  -  STr(H)  -  r(H,®)]  -  c,(H) 

xe» r  (5.3) 

(b)  V°  ■  ~{H*  -  H4®*)  •  ~  E®* 

(c)  •  o,(H)  Y  where  Y  Is  all-pass. 

Proof:  (a)  If  we  take  k«0  In  Lemma  4  then  H*o  Is  the 
optimal  zeroeth  order  Hankel  approxlmant  of  H(s)  and 

0|  (h)  -  olr(H)-  r(H4®)]  .  3H-H,°!U  . 

Since  y  is  all-pass 

oj(h)  ■  3h-h,o3..bh*-h4ou..ji  i 1  (h* - h*°*)!|.  . 

gn-HjOJl.  *  oi(H)  being  finite  H-H4®  Is  proper  and  H 
being  proper  fi4°  Is  proper.  The  order  of  Ha®  being 
zero,  Ha®  has  no  l.h.p.  poles.  Hence  H#®  has  no  r.h.p. 
poles  and  X®  £  H4°*  Is  stable  and  proper.  Thus  x°€H“ 

yields  oi(H)  -  (H*  -X®  )J*.  Me  claim  that  x® 

achieves  the  minimum  In  (5.3).  For,  If  there  existed 


Proof: 


r*(s.2.)|n  , — ] — _  (4  13)  achieves  the  minimum  In  (5.3).  For,  If  there 

‘I*  ’’  Js..,TT^7 


r"*ls.,,  ■  '('ti.  >«<**• 


Then  from  Lemna  3  (b)  Y°  is  a  solution  of  (4.12)  Iff 
it  is  a  solution  of  (4.8).  Hence  problems  3  and  4 
have  the  same  solution  YO.  „  ,  _ 


5.  Zeroeth-Order  Hankel  Norm  Approximation 

The  solution  to  the  optimization  problem  4  has 
been  found  by  Kung  and  Lin  (9.  10]  as  an  Intermediate 
Step  In  the  solution  of  the  optimal  Hankel  norm  approx¬ 
imation  problem.  With  every  transfer  function 
G(s)  €  Lm  there  is  associated  a  Hankel  operator  r(G) 
defined  as 

(T(G)u) (t)  •  Gc(t*i)  u(x)dx  o<t<«  (5.1) 

where  Gcft)  ijd'UGfs)]. ,  [-].  being  the  projection 
operator  which  retains  only  the  stable  part  of  the 
partial  fraction  expansion  of  its  argument.  The  rank 
of  r(G)  Is  the  number  of  Its  nonzero  singular  values. 

It  equals  the  order  of  [G(s)]_  which  is  defined  as  the 
dimension  of  Its  minimal  state  space  realization.  It 
equals  the  number  of  l.h.p.  poles  of  G(s).  In  the 
Hankel  norm  model  approximation  problem  the  objective 
Is  to  obtain  a  lower  order  approximation  Gets)  to  G(s) 
so  that  the  Hankel  norm  of  the  error  E(s)  t  G(s)  -  6«(s) 
is  small.  The  Hankel  norm  of  G(s)  1*  defined  as  the 
largest  singular  value  of  the  associated  Hankel  opera¬ 
tor  r(G). 

The  following  result  from  kung  and  Lin  [9,  10] 
relates  the  bounds  on  the  error  to  the  order  of  the 
approximation. 

Lemna  4:  If  [Ga(s)]„  has  order  K  m  0,  then 

(8)  *  otr(G)  -  r(G,)]  *  BG-  6,3.  (5.2) 

where  0k«l(6)  denotes  the  (k+l)st  largest  slnoular 
value  of  r(G).  Further,  there  exists  a  6a°(s)  of  order 
[G4°(s)].  ■  k  for  which  equalities  are  obtained  In 
(5.2)  and  E®  A  6  -  •  ok+,(G)  I  where  C  Is  all-pass 

and  rational.  G§®  Is  called  the  optimal  kth  order 
Hankel  approxlmant  of  G(s)  and  E®  the  corresponding 
minimal  error. 

Me  will  make  use  of  this  Lemna  to  show  that  a  so¬ 
lution  to  problem  4  can  be  obtained  via  the  optimal 
zeroeth  order  Hankel  approximation  of  H(s). 


JfH-  x*JL  <  o, (K)  and  H4  a  [X*].  is  a  zeroeth  order  ap- 
proximant  of  H(s)  which  violates  Eq.  (5.2)  of  Lemna  4. 

Thus.  o,(H)-a£  (H*-X®)3.«  Min  (H*-X);|. 

•  X€*T 

(b)  Y3  ■  £  (H*  -  X°)  -  £  (H*-  H,®«)  *  —■  E°\ 

(c)  Since  E®  •  o,(H)E  where  I  Is  all-pass  yo.il  E®*  ^ 

o,(H)Y  where  Y  £  £  1s  *)1-P***  since  V  ,nd  ^ 

are  all-pass.  Q.E.D. 

Hence  we  have  obtained 

Y°.£e®\  (5.4) 

the  solution  to  problem  4,  in  terms  of  the  minimal 
zeroeth  order  Hankel  error  E®  of  H(s).  It  gives  us 
the  optimal  weighted  sensitivity 


1®  •  Bj.*1 


Since  Br,  Dx  are  all-pass,  T®  Is  of  the  form 

T®  ■  o^(H)  T  (5.6) 

where  t  Is  all-pass.  Mote  that  all  our  optimal  solu¬ 
tions,  Ha®,  X®,  E®,  Y®.  T®  are  rational  [9,  10]  which 
justifies  the  assumption  we  made  In  defining  problem  2. 

Remark:  The  all -pass  form  of  the  optimal  weighted  sen¬ 
sitivity  Is  In  agreement  with  the  S1S0  results  obtained 
In  [3,4].  It  Is  likely  that  the  optimal  compensator 
given  by  (4.10)  will  turn  out  to  be  Improper.  In  that 
case  we  can  approximate  the  Improper  compensator  by  a 

5  roper  one  by  using  high  frequency  attenuation  as  in 
3,4]  and  the  approximation  can  be  made  to  yield  sen¬ 
sitivities  arbitrarily  close  to  the  optimal  over  any 
finite  bandwidth. 

To  compute  o,(H)  and/or  the  optimal  zeroeth  order 
Hankel  approxlmant  of  H(s)  one  may  use  the  algorithm 
given  by  Lin  and  Kung  [9,  10].  The  algorithm  involves 
computing  the  largest  singular  value  of  r(H)  and  solv¬ 
ing  a  matrix  polynomial  equation  and  an  algebraic 
rlccatl  equation.  In  order  to  obtain  h(s)  we  have  to 
compute  Br  and  V*  (Eq.  (4.4),  (4.13)).  This  requires 
left  and  right  coprime  fractional  representation  of  P, 
solving  a  left  Bezout  Identity  and  inner-outer  factori¬ 
zations  of  stable  rational  matrices.  An  algorithm  for 
obtaining  a  coprime  fractional  representation  of  ra¬ 
tional,  proper  matrices  and  solving  the  corresponding 
Bezout  Identity  Is  given  In  [13].  An  algorithm  for 


computing  inner-outer  factorizations  is  given  in  [14]. 

If  one  desires  to  compute  only  the  minimal  value 
oi(H}  of  the  cost  (4.12).  another  simple  approach  is 
afforded  by  the  identity  11] 

Oj (H)  «  (Xmax  (M^))*1 

where  xmax(-)  denotes  the  greatest  eigenvalue  and  Mo 

and  ^  are  the  observability  and  reachability  gramiians  12] 

of  any  state-space  realization  of  H(s)  [16]. 

6.  Conclusions  13] 

The  main  contribution  of  this  paper  lies  In  estab¬ 
lishing  a  link  between  the  multivariable  L"  sensitivity 
optimization  problem  and  the  optimal  Hankel  approxima-  14] 
tion  problem.  As  in  [3,4]  we  have  used  the  fractional 
representation  approach  to  parameterize  the  set  of  real¬ 
izable  sensitivity  functions.  With  this  parameteriza-  15] 

tion,  the  problem  of  minimizing  sensitivity  over  the 
set  of  stabilizing  compensators  leads  to  an  optimal  iT-  16] 
interpolation  problem. 

We  have  solved  this  problem  by  solving  the  equiva¬ 
lent  optimal  zeroeth  order  Hankel  approximation  problem. 

The  optimal  sensitivity  is  related  to  the  optimal  Hankel 
error  in  a  simple  way,  and  the  optimal  Hankel  error  can 
be  computed  by  using  an  algorithm  of  Lin  and  Kung  [9,10] 

The  approach  taken  in  [4,5,6]  to  solve  the  optimal 
iT- interpolation  problem  is  based  on  Sarason's  theory 
[15]  and  makes  use  of  the  Nevanlinna  algorithm  to  com¬ 
pute  the  optimal  sensitivity.  If  we  exploit  the  rela¬ 
tionship  established  between  optimal  iT-interpolation 
and  optimal  Hankel  approximation  we  can  use  the  Nevan¬ 
linna  algorithm  for  computing  optimal  zeroeth  order 
Hankel  approximants  as  well.  This  points  to  the  possi¬ 
bility  of  using  the  Nevanlinna  algorithm  in  higher  order 
Hankel  approximation  problems  as  well,  assuming  that 
these  problems  could  be  shown  to  be  equivalent  to  some 
iT-interpolation  problem.  This  possibility  is  not  a 
remote  one  given  the  fundamental  relationship  that 
exists  between  iT-interpolation  and  Hankel  approximation 
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